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Abstract 



We make use of the global symmetries of the Yang-Mills theory on the lattice to de- 
sign a new computational strategy for extracting glueball masses and matrix elements 
which achieves an exponential reduction of the statistical error with respect to standard 
5^ I techniques. By generalizing our previous work on the parity symmetry, the partition 

function of the theory is decomposed into a sum of path integrals each giving the con- 
tribution from multiplets of states with fixed quantum numbers associated to parity, 
charge conjugation, translations, rotations and central conjugations Z^. Ratios of path 
integrals and correlation functions can then be computed with a multi-level Monte 
Carlo integration scheme whose numerical cost, at a fixed statistical precision and at 
asymptotically large times, increases power-like with the time extent of the lattice. The 
strategy is implemented for the SU(3) Yang-Mills theory, and a full-fledged computa- 
tion of the mass and multiplicity of the lightest glueball with vacuum quantum numbers 
is carried out at a lattice spacing of 0.17 fm. 



1 Introduction 



The existence of glueballs is a distinctive property of quantized non-Abelian gauge the- 
ories [iHS]. Since the first Monte Carlo simulations of lattice field theories, glueballs 
have been the focus of many studies (see Ref. fJJ and reference therein) . The main diffi- 
culty for isolating their contribution in correlation functions, and thus computing their 
masses and matrix elements, was identified almost immediately: the signal-to-noise ra- 
tio of suitable two-point correlation functions [5j decreases exponentially with the time 
separation of the sources, and in practice it is very difficult to find a window where 
statistical and systematic errors are both under control [6l[7]. A widely used strategy 
to mitigate this problem is to reduce the contamination from excited states in the cor- 
relators by constructing interpolating fields with a small overlap over them [8])[9j. The 
lowest energy is then extracted at short time-distances by assuming a negligible con- 
tribution from excited states, sometimes also with the help of anisotropic lattices [lOj . 
The most comprehensive studies of glueball masses and matrix elements performed with 
these techniques can be found in Refs. [Sl IllffTS] . This approach is not entirely satisfac- 
tory from a conceptual and a practical point of view. The problem of the exponential 
degradation of the signal-to-noise ratio remains unsolved, and the functional form of the 
sources is usually optimized so that correlators or combinations of them show a single 
exponential decay in the short time-range allowed by the statistical accuracy. A solid 
evidence that a single state dominates, i.e. a long exponential decay over many orders 
of magnitude, is thus missing. In most of the cases the computation of more involved 
correlation functions remains inaccessible. 

Recently we proposed a new approach to solve the problem [T6lfT7] (see also [T8lfT9] ). 
By using the transfer matrix formalism, it is possible to introduce in the partition func- 
tion projectors which select the contributions from states with a given set of quantum 
numbers only. The composition rules of the corresponding transfer matrix elements can 
then be exploited to implement a hierarchical multi-level integration procedure. By it- 
erating over several levels, the numerical cost for computing ratios of partition functions 
or correlation functions grows, at asymptotically large times, with a power of the time 
extent of the lattice rather than exponentially. 

The aim of this paper is to design a new strategy for computing the lightest glueball 
masses, multiplicities and their matrix elements in each super-selection sector of the 
Hilbert space of the Yang-Mills theory. This is achieved by generalizing the analysis 
in Ref. |16] to all discrete symmetries of the lattice theory in finite volume: parity, 
charge conjugation, rotations, translations and central conjugations Zff. We assume 
some familiarity with Ref. [16], especially in Section 6, where the numerical algorithm is 
described. Being based on the transfer matrix formalism, the approach is rather general 
and can be applied to a wide class of bosonic theories. 

We implement the strategy in the SU(3) Yang-Mills theory on lattices with a 
spacing of roughly 0.17 fm, spatial volumes up to 5fm'^, and time extent up to 2fm. 
A full-fledged computation of the mass and multiplicity of the lightest glueball with 
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vacuum quantum numbers is then carried out. The algorithm behaves as expected, and 
the multi-level integration scheme achieves an exponential enhancement of the signal- 
to-noise ratio for the quantities considered. As a result we can extract with confidence 
the contribution from the lightest glueball and, from a lattice with a time extent of 
2 fm, determine its mass with a precision of few percent. Its multiplicity turns out to be 
one within statistical errors, and the latter are small enough to exclude all other values 
allowed by the underlying group theory. 



2 Preliminaries and basic notation 



For definiteness we focus on the SU(3) Yang-Mills theory discretized by the standard 
plaquette Wilson action. The theoretical discussion, however, can be applied to other 
bosonic field theories, and it is mostly independent on the discretization details. The 
theory is set up on a finite four-dimensional lattice of volume V = T x with a 
spacing a and periodic boundary condition^. The gauge-invariant action is defined as 
(all unexplained notation can be found in Ref. [16]) 



l-^ReTr{c/^^(x)} 



(2.1) 



where the trace is over the color index, and /3 = Q/g^ with being the bare coupling 
constant. The plaquette is defined as a function of the gauge links U^{x) as 

U^,{x) = U^(x) U,{x + ft) Ulix + u) Ul{x) , (2.2) 

with ^, 1/ = 0, . . . , 3, /i is the unit vector along the direction /i, and x is the space-time 
coordinate. The path integral is defined as usual 
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where VU is the invariant Haar measure on the SU(3) group, which throughout the 
paper will be always normalized so that J DU = 1. The "coordinate" basis in the 
Hilbert space of the theory is the set of vectors which diagonalize the field operators at 
all spatial points, so that on a given time-slice of the lattice 

Ufc(x)|C/,,) = [/fc(xo,x)|C/,o) . (2.4) 

As the notation suggests, the operator eigenvalues Uk{xo,x) are identified with the 
spatial links on the time-slice xq at the spatial coordinate x. The matrix elements of 
the transfer operator T are 



U,,+i\TFg\Uxo) , (2.5) 



'^Throughout the paper dimensionful quantities are always expressed in units of a unless explicitly 
specified. 
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where Pg is the projector onto the gauge invariant states. In the foUowing, for notational 
simphcity, this projector is always included in the definition of T. For the Wilson action 
the transfer matrix elements are I20fl23l 



T 



D[J7] = JJPfi(x) 



(2.6) 



where the explicit form of the Lagrangian L can be found in Appendix [A) From this 
formula one can define the transfer matrix elements associated to a given thick time- 
slice, i.e. the ensemble of points in the sub-lattice with time coordinates in a given 
interval [xo,yo] and bounded by the equal-time hyper-planes at times xq and yo; as 



«'0=2^0 + l 



D3[^]=n n^f^^w 



X fc = l 



(2.7) 

By identifying the gauge transformations in Eq. (12. 6p with the links in the temporal 
direction, the path integral can thus be written as 



T-l 

n D3[C/.o]T 

2;o=0 



which corresponds to 



Z = Tr |t^| 
with the trace being over all gauge invariant states. 



(2.^ 



(2.9) 



3 Decomposition of path integrals 

The partition function Z can be decomposed into a sum of path integrals each giving the 
contribution from multiplets of states which transform as an irreducible representation 
of a symmetry group of the theory. The phase space of the theory can indeed be 
divided into regular representations of the group, which in turn can be decomposed 
into irreducible ones by applying the standard group-theory machinery [241125] . The 
invariance of the transfer operator T under the group transformations then guarantees 
the decomposition of the partition function. An analogous decomposition applies to the 
correlation functions of the theory. 

To carry out this program in detail, let us assume that the theory is invariant under 
the transformations of a generic discret^ symmetry group Q of order g. Its elements 
TZi, with i = 1, . . . , g, act on a generic vector \U) of the coordinate basis of the Hilbert 
space as 

r{n,)\u) = |c/^o , (3.1) 

^The generalization to Lie groups is straightforward. 
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where U''^' is the gauge field obtained by applying the group transformation TZi to the 
original one (see next section). The vectors form a regular representation, i.e. 

f(7^0|f/^^) = ^r[f (7^,)|c/^0 , (3.2) 



1=1 



where 



^;f(7^, 



1 if 111 = TZiTZj 
if 7^/ / TZiTZj 



(3.3) 



The non-equivalent irreducible representations T^f^\TZi) are labeled by ^ = 1, . . . , A'r, 
where Nr is the number of classes in which the group can be divided, their dimensions 
satisfy 



Nr 

E 

,1=1 



n, 



9 , 



(3.4) 



and their characters are given by 



X 



(3.5) 



For a given realization of an irreducible representation fi, the "projector" operators are 
defined as usual by 



pM = !hiY^r\>;'>*{TZ,)r{TZ, 
^ i=i 



and they satisfy 



p(m) t _ p(m) pH _ sTAti/ s: p(m) p(A') _ ?| 



(3.6) 



(3.7) 



M=i i=i 



In particular pj^'' is the projector onto states of the Hilbert space that transform as 
the component j of an irreducible representation /i. The path integral can then be 
decomposed as 

Nr 

Z = Y, Z^^^ , Z^"^ = E {t^pJ? } • (3-8) 
By inserting a complete set of eigenstates of the Hamiltonian we can also write 



-EqT 



1 + 



^(M)=e-i?oT^^Me-^^"'^, (3.9) 



where = 1 corresponds to the invariant singlet representation. In these expressions 
Eq is the vacuum energy, £'m^ are the energies (with respect to the vacuum one) of 
eigenstates in the sector //, and are the corresponding multiplicities. The latter are 
integers and positive since for the Wilson action the transfer operator T is self-adjoint 
and strictly positive [21] . 
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3.1 Partition function 



The transfer matrix elements among states belonging to irreducible representations can 
be written as 



T 



jj 



n. 
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J]^(.^)(7^.)«VllT|f/. 



(3.10) 



For a thick time-slice the analogous ones are defined by exploiting the orthogonality of 
the projectors 



which also implies 



D3[^.o]T 



T 



T [Uy(^,Uxo] T [Uyf,,UxQ] J T [C/yg, ?7j,g_i] 



a:oJ 



-S[U] jj 



(3.11) 



(3.12) 



where the integration is restricted to the active links of the thick time-slice. The rela- 
tions (13. lip and ()3.12p are the basic building blocks for the practical implementation 
of the multi-level integration algorithm described in Section [6l Finally, by repeatedly 
using Eq. ()3.1ip . it is possible to rewrite the path integral for a given representation as 



T-l 



(m) 
jj 



(3.13) 



It is interesting to notice that even though the transfer matrix formalism inspired the 
construction, the above considerations hold independently of the existence of a positive 



self-adjoint transfer operator. The insertion of T 

,(m) 



(m) 
jj 



UxQ + l ) Uxq 



in the path integral 



plays the same role as Pjj in Eq. (j3.8p . i.e. it allows the propagation in the time 



direction of states belonging to irreducible representations /i only (component j). 



3.2 Correlators of composite operators 

We are interested in irreducible tensor operators, which under the group transform as 



(3.14) 



1=1 



that are diagonal in coordinate space, i.e. 



O'^'^^IU) = o'-''\u)\U) (3.15) 



where o\^\u) is a functional of the links. The transformation rule (I3.14P implies that 



o(.^)(c/^0 = ^^^;)*(7^,)o^^^([/). 



(3.16) 



1=1 
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By using the Wigner-Eckart theorem, we can define the reduced matrix elements cor- 
responding to an operator insertion on the time-sHce xq inside a thick time-shce of size 
d as 



(Mi)(M2) 



U, 



zo+d: 



n-- 



Ml k 
h 



(3.17) 



^1 1=1 ii=lZ2 = l 



where zq < xq < ZQ + d, ^ j are the Clebsch-Gordan coefficients, k = I, . . . , n^^^ 

and n'l^i^ is the number of times that the irreducible representation fii appears in the 
direct product representation T^^^^ (8) r*^'^^)^ Thanks to the orthogonality properties of 
the Clebsch-Gordan coefficients, a two-point correlation function can be written as 

(wq <yo<wo + d) 

j;(of)(yo)of)*(xo)) = - / D3[t/^o+jD3[[/^o]D3[t/.o+,]D3[C/.o]D3[C/o] x 
i=i 



,MM2 
'Ml 



X T 



Ml>A*2 Jl,J2 
(M2) 



(Mi)(M2) 



M2 



fc=l 



(Ml)(/i2) 



(3.18) 



J2J2 



UzQ , UzQ+d 



n(Mi) 

iiii 



where ZQ+d < wq. The generalization to n-point correlation functions is straightforward. 



4 Discrete symmetries of the Yang-Mills theory 

The SU(3) Yang-Mills theory on a lattice of finite volume with periodic boundary 
conditions is invariant under parity, charge conjugation, translations, rotations and 
central conjugations Z^. In this section we set the notation for these groups, and we 
briefly review the properties which are relevant for the paper. 



4.1 Parity 

The group is of order 2. A regular representation is two-dimensional, and it is spanned 
by the vectors 

= \U) , \U^^) = \Un , (4.1) 

where on a generic time-slice xq 

U^{xo,^) = ul{xo,-^-k) . (4.2) 

The two irreducible representations of dimension one are r(^)(7^i) = ibr(^)(7^2) = 1; 
where the phase convention is the same as in Ref. |16j . 
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4.2 Charge conjugation 

The group is of order 2. A regular representation is two-dimensional and it is spanned 

by 

= \U) , \U'^^) = \U^) , (4.3) 

where 

C/f(xo,x) = C/,*(xo,x). (4.4) 
The two irreducible representations of dimension one are r(^)(7^i) = ±r(^)(7^2) = 1- 

4.3 Translations 

The group of translations is a direct product of three Abelian groups, one for each 
space direction. Its elements are labeled by a three dimensional vector of integers 
m = (tti-i, 7712, w-a), with = 0, . . . , L — 1, where each component labels the elements 
of the Abelian group in the corresponding direction. A regular representation is L^- 
dimensional and is spanned by 

IC/^-) = If/"") (4.5) 

where 

C/r(xo,x) = t/fc(xo,x-m) . (4.6) 

Since the group is Abelian, each element forms its own class and there are inequivalent 
irreducible representations of dimension 1 

^(P)(7^^,) = e^P-™ (4.7) 
27r 

which are labeled by momentum vectors p = — [ni, n2, na], with = 0, . . . , L — 1. 

Ij 

4.4 Rotations 

The octahedral group is of order 24. Its elements are listed in Appendix iBl They form 
5 equivalence classes. A regular representation is 24-dimensional and is spanned by 

= , i = l,...,24 (4.8) 

where 



[/fc(xo,Ri^x) (Ri)/,fe>0 
f^r(^o,x) = <' . (4.9) 

I ^72(xo,R-'x- A:) (R,)ftfe<0 

The inequivalent irreducible representations are two singlets Ai and A2, one doublet 
E and two triplets Ti and T2. Their expressions and their characters are given in 
Appendix [B] while the Clebsch-Gordan coefficients can be found in Ref. 
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4.5 Central conjugations 

The presence of this symmetry, first described by 't Hooft in Ref. [26] (see also Ref. [27j 
for a review), is due to the choice of periodic boundary conditions and it disappears 
in the infinite volume limit. The group is a direct product of three Z3, one for each 
spatial direction. It is of order 27, and its elements are labeled by a three dimensional 
vector of integers v = (ui, 1^2, 1^3), with i/j = 0,1,2, where each component labels the 
elements of the Abelian group in the corresponding direction. A regular representation 
is 27-dimensional and is spanned by 

pTlu^ = \u'') (4.10) 

with 

C/fe"(xo,x) = A,(x)C/fc(xo,x)At(x + A:) , (4.11) 

and 

A^(x) = e*ir('^-^)^ , A^(x + LA;) = e^^'^'=A^(x) (4.12) 

where is a 3 x 3 diagonal matrix with elements Waa = (1 — 3 ^as)- Since the group is 
Abelian, each element forms its own class, and there are 27 non-equivalent irreducible 
ones of dimension 1 

^('=)(7^^) = e^''"^ (4.13) 
27r 

which are labeled by the electric flux vectors e = — [ei, 62, 63], with a = 0, 1, 2. 



5 Glueball masses, multiplicities and matrix elements 

Dynamical properties of glueballs can be extracted from ratios of partition functions and 
correlators, which in turn can be computed efficiently by judiciously putting together the 
tools developed in the previous sections. As usual, the various super-selection sectors of 
the Hilbert space are identified by the quantum numbers associated to a complete set 
of operators which commute among themselves and with the transfer operator. In the 
zero- momentum sector (p = 0), glueball states are classified by their transformation 
properties under rotations (/^, j), parity (V) and charge conjugation (C). Moreover in 
finite volume a null electric fiux vector (e = 0) identifies the "physical sector" of the 
theory, i.e. the one which survives in the infinite volume limit. The corresponding 
projectors are given by 

p(M,-p,C) ^ p(M)p{7:')p(c)p{p=0)p(e=0) ^ ^^^^-^ 

and the partition functions read 

z(^'^'^) = ^Tr{t^pg'^''^} . (5.2) 
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The transfer matrix elements associated to a given thick-time shce are defined as 

(5.3) 



rjn(M.T'.C) 



33 



-^i:R.i.D.c 



2592 ^ ^j;R».p.c 



where 



and for each group the sum is on ah its elements. The corresponding path integrals are 
finally given by 

^(.,.,c) = ^ / J] D3[C/.J Tg'-'^) [Ux,^,, Ux,] . (5.5) 

j=l-^ xo=0 

In each sector, thanks to Eqs. ()3.9p and with the exception of the vacuum one (see 
below), the contribution of the lightest glueball is the leading exponential in T. Its mass 
and multiplicity can thus be extracted from the ratio Z^^'''^''^^ /Z^^^'^'^^ computed at 
large enough values of T, where the contamination from heavier states can be neglected. 

The sector (//,P,C) = (Ai,+,+) is special because the vacuum contribution dom- 
inates the partition function at large T. In this case the mass and the multiplicity of 
the lightest glueball can be determined by a two step procedur^. First, one defines the 
projector onto physical states with non-zero momenta and positive charge conjugation 

p(p,+) ^ p(C=+)p(p)p(e=0) ^ ^^ Q-^ 

where, for instance, p = [27r/L, 0, 0]. The vacuum does not contribute to the correspond- 
ing partition function Z^P'"*"-* which, at large T, is dominated by the lightest glueball 
state even under charge conjugation. Its energy and multiplicity can be extracted from 
the large T behaviour of Z(P'+) The mass can then be determined, up to 0{a'^) 

effects, by using the continuum dispersion relation. In the second step one defines the 
projector 

p{Ai,+,+) ^ p(p=0)p(e=0) jji^ _ p(Ai,+,+)| ^ ^^^rj^ 

computes z\^^'~^'~^'^ defined analogously as in Eq. (|5.5p . and extract the mass and the 



multiplicity of the lightest state. If, in the continuum limit, the mass turns out to be 
heavier than the one computed in the first step, the latter is the mass of the lightest 
glueball with vacuum quantum numbers, and it is also the mass gap of the theory. 

By generalizing Eq. (jS.lSp to the case of multiple quantum numbers, correlation 
functions of composite operators can be calculated analogously to ratios of partition 
functions discussed above. It is interesting to notice that the approach described here 
allows to determine efficiently also the mass and the multiplicity of states in the "un- 
physical" sectors (e / 0) [28j . 



^This is the simplest example of how the procedure described in this paper can be generalized to the 
computation of excited states in each symmetry sector. 
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Figure 1: Schematic description of the partitioning of the lattice in the multi- level 
integration algorithm. 

6 Numerical algorithm 

The formula in Eqs. ()5.3p and (15. 5p require some manipulation before they can be 
implemented in a numerical simulation. 

6.1 Thick-time slice with fixed quantum numbers 

In the computation of the transfer matrix elements (normalized to the standard one) 
given in Eq. (|5.3p . the basic building blocks are the ratios 



T 






T 


UyQ , UxQ 





Once written as in Eq. (j3.12p . they are computed numerically through the telescopic 
algorithm described in Section 4.1 of Ref. [16|, the latter being generalized to all possible 
combinations of group transformations. Their calculation is the most expensive part of 
the multi-level procedure. It is therefore worthwhile to optimize on the number of them 
and/or on their numerical precision. To this aim it is relevant to notice that a weighted 
average appears in Eq. (jS.Sp . with the weights given by the appropriate products of 
characters. The larger is the number of significant terms averaged over, the lower is the 
statistical precision required on each term so to achieve a given accuracy on the sum. 
In the case of the singlet under all symmetry groups, for instance, the precision on each 
thick-time slice ratio can be reduced proportionally to the square root of the number 
of addenda. A large number of terms therefore does not automatically implies a more 
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expensive numerical computation. The statistical error on each ratio ()6.ip . however, 
cannot be arbitrarily large for these arguments to apply. Therefore it also pays off to 
reduce the number of ratios to be determined. This can be achieved by promoting 
each group index to a stochastic variable with integer values. To each of the index 
combinations is thus associated a probability distribution IlRi^p^c,m,u so that 



Hi ,p,c,m,i> 



1 . 



Moreover if we sum over a subset of indices, for instance Rj, then 



n 



p,c, m,i/ 



En 



(6.2) 



(6.3) 



Ri 



is still a probability distribution in the remaining variables. If, for instance, we choose 
^B.i,p,c,m.,u = const, then the thick-time slice with fixed quantum numbers can be written 
as 



33 




T 







Ri ,p,c,m,i> 



T 



TTRi,p,c,m,i/ 



,U 



XO 



Uyo , Uxo 



(6.4) 



The sum on the r.h.s can be computed stochastically by extracting, for each thick-time 
slice, a series of "configurations" for the set of indices (Rj,p, c, m, iv) and averaging 
over them. There is clearly a great freedom in choosing the best sampling procedure. 
An efficient one is to project exactly on the quantum numbers which eliminate the 
contributions from states lighter than the one of interest, while treating stochastically 
the indices which project out heavier states. One can also extract the indices with 
a distribution different from a constant one, or implement a more involved sampling 
procedure. Once the estimate of the l.h.s. of Eq. (|6.4p is inserted in the multi-level 
algorithm (see below), the final result is independent from the particular procedure 
implemented and from the statistical accuracy on the thick time-slice matrix elements. 
The algorithm is design to be always exact, but the variance and therefore its efficiency 
depends on details of the implementation. 



6.2 Multi-level integration scheme 

The ratio of partition functions Z'^^''^'^^ jZ can then be calculated by implementing the 
hierarchical two-level integration formula (see Figure [T]) 



where 



S{U\ pCM.-P.c) 



p(M.-P,c) 

j;m,d 



m-l rp(At,-P,c)r7-T TT 1 

■j-j- ^jj ['-^xo+{i+l)-d,'-^xo+i-d\ 

i=0 



xo+ii+l)-d7 Uxo+i-d 



(6.5) 



(6.6) 
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with m > 1, uq = (xq + m ■ d). In each ratio on the r.h.s of Eq. (j6.6p the stochastic 
indices are generated, independently on each thick-time sHce, with chosen distribution 
nRi,p,c,m,i/- The procedure can be generahzed easily. For a three-level scheme, for 
instance, each ratio can be computed with a two-level scheme by exploiting the compo- 
sition rule in Eq. (|3.1ip . While the result does not depend on the particular integration 
scheme implemented, its statistical error does. The algorithm therefore requires an 
optimization which allows one to exploit the expected spectral properties of the the- 
ory. For instance if d is chosen large enough, i.e. larger than 1/Tc with Tc being the 
critical temperature, only a few of the physical states give a sizeable contribution to 
each ratio T^^'^''^\Uxg+d,Uxo]/T[Uxg+d,UxQ]- The latter is therefore expected to be of 

order e i , the magnitude of the product is of the order of e i for each 

configuration of the boundary fields, and the statistical fiuctuations are reduced to this 
level ^EH- 

6.3 Correlators of composite operators 

The formula in Eq. (I3.18[) and its generalization to the case of multiple quantum numbers 
can be implemented by inserting the sources in a thick time-slice and then following 
the same steps that lead to Eq. (j6.5p . It is important to stress, however, that once 
all the thick time-slice ratios T^-^' ' \UxQ+d^Uxo\/T^[Ux(i+diUxo] have been computed, 
only the first integral of the telescopic expansion in Eq. (4.2) of Ref. [16] needs to be 
re-done for the thick time-slices where a source is inserted. The extra numerical burden 
is therefore negligible. As in the previous subsection, the statistical error associated to 
the estimate of the correlation function is comparable to the signal if the algorithm is 
properly optimized. 

7 Numerical results 

To demonstrate the practical feasibility of the strategy proposed in this paper, we 
have carried out a full-fiedged computation of the mass and multiplicity of the lightest 
glueball with vacuum quantum numbers. We have simulated the SU(3) gauge theory 
at /? = 6/^0 = ^•'''j which corresponds to a spacing of 0.17 fm if the reference scale 
ro = 0.5 fm is used to calibrate the lattice [29]. The spatial lengths of the lattices 
are 1.4 and 1.7 fm, while their time dimension extends up to 2 fm. A list of the runs, 
the number of configurations generated and some details of the multi-level algorithm 
implemented is reported in Table [H 

The primary quantity that we have calculated numerically is the ratio of partition 
functions Z^P'^^^/Z defined in section [5l with p2,3 = 0, pi = (27r/L)ni and ni = 1,2. 
The thick time-slice transfer matrix elements associated to the projector in Eq. (jS.Gh 
have been computed as described in section [6j the sum on the group indices of charge 
conjugation and translations along direction 1 has been done exactly, while the one 
over the remaining indices (p2)P3,*^) has been carried out stochastically by extracting 
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Lattice 


L 


T 


A'conf 




d 


Ai 


8 


4 


50 


2 


4 


A2 




5 


50 


2 


5 


A3 




6 


100 


2 


3 


A4 




8 


100 


2 


4 


As 




12 


50 


3 


{3,6} 


B3 


10 


6 


50 


2 


3 



Table 1: Simulation parameters: A'^conf is the number of configurations of the uppermost 
level, A^i ev is the number of levels and d is the thickness of the thick time-slice used for 
the various levels. 



between 9 and 64 random "configurations" of indices with a flat probability. The results 
are collected in Table [21 and those of the A lattices are plotted in Figure [2] as a function 
of T (left panel). The data show a clear exponential decay of the ratio Z^^'^^ / Z^^'~^^ 
over more than 6 orders of magnitude. 

We fit the results of the A series to a single exponential, i.e. 



In 



A-BT , 



E 



(p,+) 



(7.1) 



This 



where A = Inw^, with ui^ being the multiplicity of the state, and B — j^^q 
function fits well the last four points, and the best fit give^^ = —0.6(4) and B = 1.15(6) 
(x^/dof = 1.5). Group theory predicts the multiplicity to be an integer between and 
3. The null value is excluded by the data, the multiplicity 1 is within 1.5(T of the central 
value given by the fit, while 2 and 3 are 3.2 a and 4.2cj away. As a further check we 
also fix the multiplicity to one of the possible integers, and we obtain x^/dof = 1.7, 4.4 
and 6.8 for = 1, 2 and 3 respectively. We therefore conclude that the data strongly 
prefer multiplicity 1, the value expected for a singlet under the octahedral group. We 
stress once again that the computation of the multiplicity is new because this quantity 
is not accessible to the standard technique. 

By imposing the multiplicity to be 1, the energy of the state can be determined at 
each T as 



E. 



(p,+) 

off 



4 in 

T 



(7.2) 



which yields the results given in Table [2] and shown in the right panel of Figure [2j It 
is interesting to notice that a precision of a few percent is reached with only 50-100 



*The small value observed at T = 6 is compatible with a 2<j statistical fluctuation, and it is responsible 
for the value of x^/dof a bit larger than expected. 
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T J.J.- 7(0, + ) 1 rv 

ijattice Zj^ ' ' / Zj 


r, 7(P, + ) 17 7(P,+) /7(0,+) i7i(P. + ) 

n\ Zj^^^ ' 1 Zj z^^' ' / z^ ' ' -^cff 


Ai 0.158(14) 

r^2 U.o ( ^\0) 

A3 0.96(4) 
A4 0.999(9) 
A5 0.963(13) 


1 1.6(3) • 10-3 1.04(21) • 10-2 1.14(5) 

1 4.5(7) • 10-^ 4.7(7) • 10-^ 1.277(25) 
1 6.6(12) • 10-^ 6.6(12) • 10-^ 1.203(22) 
1 4.1(16) • lO"'^ 4.3(17) • lO"'^ 1.22(3) 


B3 1.03(4) 


1 1.0(3) • 10-3 1.0(3) • 10-3 1.15(5) 

2 0.94(25) • 10-^ 0.92(25) • IQ-'' 1.55(5) 



Table 2: Results for ratios of partition functions with momenta p = [27rni/L, 0, 0]. The 
effective energy E^^^^ is defined as in Eq. ()7.2p . 

configurations of the uppermost algorithmic level. To be on the conservative side, we 
take as our best estimate for the energy the value at T = 12 reported in Table [2j 
The result of the fit to a constant of the last four data points gives 1.233(14), a value 
compatible with our best one but with half the statistical error. 

Finite volume effects in the energy values are expected to be exponentially sup- 
pressed at asymptotically large values of L. Lattice B3 serves the purpose of assessing 
their magnitude. It has the same lattice spacing of the A series but a linear extension of 
L = 10. The results for ni = 1, 2 are reported in Tabled and are plotted as a function 
of the momentum squared in Figure [3j The dashed line is a linear interpolation of the 
two black points (circles) of the lattice B3, while the red point (square) is our best result 
for the A series. It is rather clear that, within our statistical precision, finite volume 
effects are not visible in our data. It is also interesting to notice that, even if the values 
of the momenta are rather large, the continuum dispersion relation is well reproduced 
within our statistical errors. 

The glueball mass can finally be extracted, up to 0{a^) discretization errors, by 
using the continuum dispersion relation 

M+ = ^(i?S'"'^)2-p2 , (7.3) 

which, in units of the lattice spacing, yields to 

M+ = 0.935 ± 0.042 . (7.4) 

This is one of the main numerical results of this paper. When converted in physical 
units by using tq, it gives M"*" = 1.08(5) GeV. We remark that the value in Eq. (I7.4p is 
fully compatible with 0.955(15), the mass computed with the standard method by using 
the same action at the same lattice spacing [12] . On the other hand its value in physical 
units is quite smaller than the continuum extrapolated one estimated, for instance, in 
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Figure 2: Left panel: ratio of partition functions Z^^'^^ /Z^^'^^ with momenta p = 
[2tt/L, 0, 0] for the A lattices. Right panel: the corresponding effective energy as defined 
in Eq. (|7.2p . The band is our best estimate, i.e. the one extracted from the lattice with 
the longest time-extension. 



Ref. |14) . This can be easily explained by discretization effects which still contribute to 
the value in Eq. (|7.4p [12]. Removing them goes beyond the scope of this paper. 

To identify the state of mass with the 0"*"+ glueball we still need to verify 
that the lightest state in the sector orthogonal to the vacuum, i.e. the one selected by 
the projector in Eq. (|5.7p . is heavier. Since the multiplicity turns out to be one, we 
can further restrict the orthogonal sector by projecting onto the singlet representations 
of the octahedral group Ai and A2. The ratio of such a partition function over the 
standard one has been computed in a dedicated run at T = 6 and L = 8, and our best 
estimate turns out to be 2.6(26) • 10""^. This result, even if compatible with zero, puts 
an upper bound on this quantity much smaller than the one expected if a state lighter 
than M+ were present. This suggests to identify with the mass of the ^4^"'' glueball. 
More simulations are needed to corroborate this result, and to have a solid numerical 
evidence of the existence of a mass gap in the theory. 
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Figure 3: The effective energy squared from tlie (red square) and tfie (black 
circles) lattices. The line is a linear interpolation of the black points (circles). 



8 Conclusions 

The relative contributions to the partition function, due to states carrying a given set 
of quantum numbers associated with the exact symmetries of a field theory, can be 
expressed by ratios of path integrals with different boundary conditions in the time 
direction. From a theoretical point of view these are very clean quantities, which have 
a finite and universal continuum limit once the bare parameters in the action have 
been renormalized. Prom an algorithmic point of view, the composition properties 
of the projectors can be exploited to implement a hierarchical multi-level integration 
procedure which solves the problem of the exponential (in time) degradation of the 
signal-to-noise ratio. 

The numerical study presented in this paper demonstrates that glueball masses 
can be extracted from these observables (for the SU(3) Yang-Mills theory) with the 
present generation of computers. We have been able to follow the exponential decay of 
a ratio of partition functions dominated by the lightest glueball with vacuum quantum 
numbers for more than 6 orders of magnitude. A fit to a single exponential in the time 
range 0.85 - 2 fm has allowed us to determine unambiguously the multiplicity of the 
lightest state for the first time. The mass has then been extracted with a precision of 
few percent at a time distance of 2 fm, where the contamination from excited states is 
negligible. 

The ideas presented here can also lead to a solid quantitative evidence of the ex- 
istence of a mass gap in the theory. To this aim a continuum limit extrapolation of 
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the numerical results is mandatory. At present, however, the scaling of the algorithm 
with the square of the number of spatial points, together with the limited numerical 
resources at our disposal, prevents us to simulate the theory at small enough lattice 
spacings for such a limit to be reliably taken. Another intriguing application is the 
computation of the thermodynamic potentials at non-zero temperature as suggested 
in Ref. [30j. No vacuum subtraction or renormalization constant is required, and the 
method can be applied at arbitrary high and low temperatures. The ratios introduced 
here may also turn out to be useful in the study of QCD-like theories, where isolating 
the contributions from single states in correlation functions may not be easy. 
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A Lagrangian of the Wilson action 



The Lagrangian L which enters the definition of the transfer matrix elements in Eq. (|2.6p 
is given by 

1 r T 1 r T 

(A.l) 



















L 




= K 













where the kinetic and the potential contributions are defined as 



K 



and 



respectively. 



W 



Uxo 



l-^ReTr{c/fc(xo + l,x)C/2(xo,x)} 



X k,l 



(A.2) 



(A.3) 



B Irreducible representations of the octahedral group 

In this appendix we report our conventions for the octahedral group and its irreducible 
representations. The group has 24 elements Rj, 3x3 orthogonal matrices (R^Ri = 1), 
grouped in five conjugacy classes. 
First class (E): 

Ri = I 1 I . (B.l) 
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Second class (C3): 

0-10 

R2 = I -1 I , R3 = I -1 I , R4 = I 1 

-10 

1 

R5 = I 1 I , Re = I 1 I , R7 = I -1 

0-10 



Rg = -1 , Rg = 1 . (B.2) 



Third class (C2): 

100\ /-100\ /-I 00 

Rio = I -1 , Rn = 1 , R12 = -10 

00-1/ \00-l/ VOOl 



n 

u 


1 

J. 


u 


n 

u 





1 

— 1 


-1 

— 1 











1 











1 


1 














-1 


-1 











1 







u 


_1 

— J. 





n 
w 


n 
u 


1 

— 1 


1 














-1 


1 











-1 











1 


1 











1 






(B.3) 



Fourth class (C4): 



100\ /00-l\ /OlO 

Ri3 = I 1 , Ri4 = 1 , Ri5 = -1 I (B.4) 

0-10/ \ioo/ \001 

100\ /001\ /O-IO 

R16 = I -1 , Ri7 = 10, R18 = 1 

010/ \-100/ \ooi 

Fifth class (C^): 

010\ /0-10\ /OOl 

Ri9 = I 1 , R20 = -1 , R21 = -1 I (B.5) 

00-1/ \00-l/ \100 

00 -1\ /-100\ /-I 00 

R22 = I -1 , R23 = 1, R24 = -1 

-100/ \oio/ \0-10 

The group has 5 inequivalent irreducible representations. 
Singlet Ai: 

r{Ri) = 1 i = 1, . . . , 24 . (B.6) 
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hi 


8G3 


3G2 




6C2 


Ai 
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1 


1 


1 


1 


1 


A2 


2 


1 


1 


1 


-1 


-1 


E 


3 


2 


-1 


2 








Ti 


4 


3 





-1 


1 


-1 


T2 


5 


3 





-1 


-1 


1 



Table 3: Table of characters for the octahedral group 



r(Ri) = i i = i,...,i2 

r(Ri) = -1 i = 13,...,24. 

r(Ri) = r(Rio) = r(Rn) = r(Ri2) = i , 

r(R2) = r(R3) = r(R4) = r(R5) = - cos (^) l - ^ sin (^) (72 

r(R6) = r(R7) = r(R8) = r(R9) = - cos (^) i + ^ (^) 

r(Ri3) = r(Ri6) = r(R23) = r(R24) = cos (|) - sin (^) cTi 
r(Ri4) = r(Ri7) = r(R2i) = r(R22) = cos + sin cji 

r(Ri5) = r(Ri8) = r(Ri9) = r(R2o) = -a^ , 

where 1 is the 2x2 identity matrix and Oi are the Pauli matrices. 
Defining triplet T\: 

r(Ri) = Ri z = l,...,24. 

Triplet Tg: 

r(Ri) = Ri i = l,...,12 

r(Ri) = -Ri i = 13,...,24. 

The characters of the various representations are reported in Table [31 



Singlet A2: 



Doublet E: 
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